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1 . 0 INTRODUCTION 


The NASA Lewis Research Center employs a general computer program (Ref- 
erence 1) for calculating the thermodynamic performance of jet propulsion 
engines. To calculate off-dusign engine performance, the user must input 
component maps. These maps define the characteristics of the various compo- 
nents ovei their full range of operating conditions, 

For advanced propulsion systems these characteristics are not generally 
known. Furthermore, the typical user of the program is not sufficiently 
knowledgable and/or cannot afford the time to do an extensive design analysis 
of the component in question. Instead he usually scales some available map. 

The objective of the study is an improved method of representing the 
turbine component when performing calculations of off-design performance for 
advanced air-breathing jet engines. This method, wliich is a computer pro- 
gram called PART, is compatible in both form and format with the cycle pro- 
gram of Reference 1 and the example map representation of Reference 2. 

Because this report contains a description of the input-output data, 
values of typical inputs, and sample cases, it is suitable as a user's manual. 
A brief description of the engineering analysis used to generate the program 
is given near the end of the report. 

The program upes turbine design point data as input to generate off- 
design values of tun ine flow-function and total-Cn-total efficiency over a 
range of pressure ratios and speeds specified by the user. A user-specified 
option will also permit calculating design point cooling flows and the corre- 
sponding change in turbine efficiency. The cooling flow subroutine, developed 
at the Lewis Research Center, is described in Reference 3. 

The Aircraft. Engine Group uf the General Electric Company has a turbine 
data base consisting of 25 turbines having design point turbine flow func- 
tions ranging from about 14 to 290. The number of stages for each of these 
turbines together with the approximate design point values of specific work 
output divided by inlet total temperature (DHQTD) and flow function (TFFD) 

'’.r;* summarized in Table I. The last seven turbines shown in the table are 
variable-geometry turbines. Table II shows the set of first stage nozzle 
area ratios for each of the seven variable-geometry turbines. Five of these 
variable g^'ornetry turbines were generated by turbine design and off-design 
computer programs similar, if not identical, to that described in Reference 4. 
Two of the turbines shown in Table II were generated from air turbine test 
carried out by the Lewis Research Center. The results of these tests are 
given in References 5 to 8. The Table II desigi.ation of the NASA test tur- 
bines have been given a trailing X. The analytical prediction of the perfor- 
mance of these turbines obtained from Reference 9 has been a trailing P. 


1 


Table I. Summary of Turbine Design Point Data. 


No. 

Data Base Name* 

No. of Stages 

DHQTD 

TPFD 

1 

HPTl-1 

1 

0.0596 

14.6 

2 

HPT 1-2 

1 

0.0705 

16.4 

3 

HPTl-3 

1 

0.0335 

88.5 

4 

HPT2-4 

2 

0.0670 

17.3 

5 

HPT2-5 

2 

0.0787 

32.4 

6 

HPT3-6 

3 

0.0810 

45.5 

7 

LPTl-1 

1 

0.0220 

45.0 

8 

LPTl-2 

1 

0.0425 

45.0 

9 

LPT2-3 

2 

0.0571 

58.5 

10 

LPT2-4 

2 

0.065 

60.4 

11 

LPT4-5 

4 

0.0665 

106,0 

12 

LPT4-6 

4 

0.0709 

134.4 

13 

LPT6-7 

6 

0.0814 

104.9 

14 

PT3-1 

3 

0.0800 

210.0 

15 

AT3-1 

3 

0.0590 

- — 

16 

AT3-2 

3 

0.0785 


17 

\T3~3 

3 

0.0635 

43.16 

18 

AT4-4 

4 

0.0499 

0 o o c 
JU t O.J 

19 

VATl-1 

1 

0.044 

99.0 

20 

VATl-2 

1 

0.060 

60.0 

21 

VAT 1-3 

1 

0.0238 

290.0 

22 

VAT1-4P 

1 

0.0328 

61.8 

23 

VAT1-4X 

1 

0.0328 

61.8 

24 

VAT2-5P 

2 

0.0636 

61.8 

25 

VAT2-5X 

2 

0.0636 

61.8 

*HPT - 

High Pressure Turbine 



LPT ■ 

- Low Pressure Turbine 



AT - 

Air Turbine Test Rig 




VAT - 

Variable Area Turbine 
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Table II. Suroraaty of Variable Geometry Turbines 
Included on Data Base. 


Turbine 

No. 

Designation 

No. of 
Stages 

First Stage ’<ozzle Area Ratios 
% 

1 

VATl-1 

1 

50.0 62.5 75.0 87.5 100.0 

2 

VAT1«2 

1 

71.0 86.0 100.0 109.0 120.0 

3 

VAT 1-3 

1 

76.0 84.0 92.0 100.0 108.0 116.0 

4 

VAT1-4P 

1 

70.0 100.0 130.0 

5 

VAT1-4X 

1 

70.0 100.0 130.0 

6 

VAT2-5P 

2 

70.0 100.0 130.0 

7 

VAT2-5X 

2 

100.0 
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tOGRAM STRUCTURE 


A flow chart' of control in the NASA parametric turbine 
program is shown :er the input has been read and processed, 
the program carriitch line analysis starting with the last 
stage of the turbs starts at the exit of the turbine stage 
(in order to avoi calculates the bucket and nozzle flow 
angles. This st«»en used to generate the stage flow and loss 
characteristics tally based correlations developed during 
the program. Suce then calculated until the first stage is 
reached. The firristics are then generated, and the stages 
stacked for each »t stage turbine nozzle area specified. If 
the turbine is ccocedure given in Reference 3 is used to 
calculate both tHquirements , and the cooled turbine effi- 
ciency. Finally j’ocessed to obtain a turbine map representa- 
tion compatible wk of Reference 1. 





Figure 1. Flow Chart Showing F3ow o 
Turbine Program. 
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Control in Parametric 





















3.0 PROGRAM INPUTS 


All of the PART inputs are of the free- fie Id format (NAMELIST) type, and 
begin in column two. Tltere ia no specified order to the inputs. The program 
initially lists the contents of the NAMELIST INPUT together with Che default 
settings of all the input variables. The user may then change as many of the 
inputs as desired. The program then echoes the updated NAMELIST. If none of 
Che inputs are changed, the program will execute the first example case and 
the user can inspect the output. The input variables together with their de- 
fault settings are summarized in Table III. 

Tlie first six input variables in Table III are used to control the number 
and values of speed, pressure ratio, and nozzle area ratios (transformed into 
nozzle angle) to be written on the output files. For example, in the input to 
the first example case shown in Figure 2, all of the corrected speed and pres- 
sure ratio arrays are used by the program, but only the first three positions 
in Che area ratio array. Note that speeds and pressure ratios are entered in 
increasing value, but that area ratios are entered in decreasing value (this 
is 30 that the nozzle angles will be written in increasing order on the output 
file). Speeds less than 10% should not be used. The input to the second ex- 
ample case shown in Figure 3 illustrates the use of the first six variables 
to limit Che size of the output files. 

Some of the design point inputs will be calculated internally by the pro- 
gram, if the user inputs the correct value to trigger Che calculation. This 
subset of inputs together with the required settings are summarized in Table 
IV. An input value equal to or less Chan zero will trigger all of the calcu- 
lations with the exception of exit swirl angle, here a value greater than 90 
degrees must be input (180 degrees is recommended). 

A minimum set of design point input would consist of NSTG and DHQTD. 

The values of TFFD and XNFTD could be inpu; as 100,0 and the resulting values 
of TFFO interpreted as percent. The user could then use the settings In Table 
IV to trigger program calculations of the remaining design point information. 
The use of the program to calculate the number of stages will frequently re- 
sult in a single-stage turbine, since the only upper limit on turbine radius 
is the limiting value of rim speed. This is not usually sufficiently restric- 
tive to require the use of additional s..ages. 

If the user wishes the program to calculate Che value of design point 
cooling flow, and Che corresponding decrease in turbine efficiency, the JCOOL 
switch in the NAMELIST INPUT should be set to 1. The program will then list 
the contents of the NAMELIST INPUTl together with the default settings of all 
the variables. The user may then change as many of the inputs as desired. 
Since the default settings of the NAMELIST INPUT are for an uncooled turbine, 
these inputs (namely, TTIN, PTIN) must be changed in order to successfully 
calculate cooling flows. The input variables for NAMELIST INPUTl together 
with their default settings are summarized in Table V. The input to Che 
second example case shown in Figure 3 illustrates the proper format for a 
cooled three-stage turbine. 
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EXEC NP'T 

ENTRY lA) El4C6£A.Fi E15.DATA DELETED 
£NTRY lA) E1426£A,(‘..^E16.DATA DELETED 


NAMELIST 

INPUT 




NSPDS » 
APCNC (1)» 

ISt 




1 


20.000000* 

30.000000* 

40.000000! 

5 

50.100000* 

60.000000* 

70.000000* 

80.000000 

9 

90,000000* 

100.000000* , 

110.000000* 

120.000000 

13 

130.000000* 

140.000000* 

150.000000* 


NPR * 

APR (I ) • 

20* 




1 

1.100000* 

1.200000* 

1 . 400000 * 

1 . 60000 ;■ 

5 

1.700000* 

1.800000* 

2.P000P0* 

2.200i.''V 

9 

£.400000* 

2.600000* 

£ . &00OO(';^ * 

'i'.Ot'P.U'. 

13 

3. £00000* 

3.400000* 

3.600000* 

3 . 1, 

17 

4 . 000000 * 

4.200000* 

4.400000* 

4. 6 Of''' 

NAR » 

ARN (!)• 

3* 




1 

1.300000* 

1.000000* 

0.700000* 

0.0 

5 - 

0.0 , 1 

,0.0 



NSTC 0 

1* 

JCOOL » 

0* 


DHOTD * 

0.032780* 

ETATTD- 

0.923000* 


TFFD * 

6t. 979996* 

XNRTD « 

193.520004* 


PSID « 

0.851100* 

XMZYn ■ 

0.373000* 


AN&SWX« 

15.200000* 

TTiiJ • 

518.669922, 


PTIN » 

14.696000* 

FARGO » 

0.0 * 


END "namelist 

INPUT 




ENTER NAMELIS 

T INPUT •? 




5 INPUT £ END 




NAMELIST 

INPUT 




NSPDS « 
APCNC (1)» 

15* 




1 

10.000000* 

20.000000* 

30.000000* 

40.00000,;' 

5 

50.000000* 

60.000000* 

70.000000* 

80.0000PO 

9 

90.000000* 

100.000000* 

110.000000* 

120.00001.1. 

13 

130.000000* 

140.000000* 

150.000000* 


NPR * 

APR m» 

Z0I 




1 

1.100000* 

1.200000* 

1.400000* 

1 .6000130 

5 

1.700000* 

1.800000* 

2.000000* 

2.200000 

9 

E. 400000* 

2.600000* 

2:800000* 

3.00000C 

13 

3. £00000* 

3.400000* 

3 . 600000 * 

3.800000 

17 

4.000000* 

4.200000* 

4.400000* 

4.600000 

NAR « 

ARN (I)« 

3* 




1 

1 .300000* 

1 . 000000 * 

0.700000, 

0.0 

5 

0.0 * 

0.0 * 



NSTC « 

1. 

JCOOL » 

0* 


DHQTD « 

0.032780* 

ETATTD» 

0,923000* 


TFFD » 

62.979996* 

XNRTD * 

193.520004, 


PSID * 

0.851100* 

XM2XD « 

0.373000* 


AN&SWXs 

15.199997* 

TTIN * 

518.669922* 


PTIN « 

14.696000* 

FARGO » 

0.0 * 


END NAMELIST 

INPUT 





NASA OUTPUT ON 1FC»15.16 
READY 


Figure 2. Input to First Example Case. 
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READY 
EXEC NPT 

JENTRY.JAJuXJ.J2i>2A .f ILEIS^IA JIELETED 
ENTRY (A) E14262A.FILE16.DATA DELETED 
NAMELIST , INPUT 


NSPDS * 





APCNC (D* 





1 

10. 000000 r 

20. 000000 r 

30.000000. 

40.000000. 

5 

^0.000000» 

_ _ 60.000000. 

70.000000. 

80.000000. 

9 

90.000000f 

100. 000000 r 

110.000000. 

120.000000. 

13 

130. 000000 > 

140. 000000 > 

150.000000. 


NPR * 

_..20t .. 




APR (D* 





1 

1. 1000001 

1. 200000 f 

1.400000. 

1.600000. 

„J5 

1. 700000 

. 1^300000 » 

2.000000. 

2.200000. 

9 ^ 

2. 400000 » 

2. 600000 > 

2.800000. 

3.000000. 

13 

3. 200000 r 

3.400000. 

3.600000. 

3.800000. 

17 

4. 00.0000 1-„ 

4.2000001 

. 4^^0000 £ . .. 

600000. 

NAR * 

3i 




ARN ( I ) * 





1 

..X*500000»> 

1 • 000000 f 

_ 0.700000. 

0.0 . 

5 

0.0 » 

0.0 . 



NSTG « 

1 r 

JCOOL * 

0. 


riWQTn s 

0.032780» 

ETATTDs 

0.923000. 


TFFD * 

62. 979996 r 

XNRTD « 

193.520004. 


PSID » 

0.851100> 

XMZXD ■ 

0.373000. 


ANGSWX* 

15.200000f 

TTIN » 

518.669922. 


PTIN * 

14.696000. 

FARGD * 

0.0 . 


END NAMELIST 

INPUT 





ENTER NAMELIST INPUT ? 

INPUT . , 

NSPDSs6 

APCNC« 20 . f 40 . » 60 . f 80 . » 1 00 . » 1 20 . 

„ NPR=5 ^ - ....... 

APRs2.0r2.5»3.0f3.5»3.8 

NAR=l 

ARN:=1.0 „ _ 

NSTG=3 
JCOOL= 1 

DHQTD= 0.0635 „„ 

ETATTD=.886 

TFFD=58.53 

XNRTD=40. 1 0 ’ 

PSID=i.5 

XM2XD=0.41 

_ftNPSWX=2.,.? _..._ 

TTIN=2500. 

PTIN=59.8' 

_FARGD=0.02 

fcEND 

Figure 3. Input to Second Example Case. 
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NAMELIST INPUTl 

_KINDOF«| 66400^0 0 r ' 

fc " « 7007000000 f FARCir^ 070 r 

YEAR s , 1980. » YEARB ■ 1980. » 

ELIFE s 0.1000000E-»-05 f _____ 

END NAMELIST "lNPUTr“ 


ENTER NAMELIST IN PUTl ? 
"CitsiPUTi 
KINDOFs864000 
SiEND 

NAMEL IST INPUT 
NSPDS * 6f 

APCNC (I)= 


1 20.000000r 40.000000> 60.000000r 

5 100.0000001 120.000000> 70.000000r 

9 90 . 000000 > 100. 000000 > L10 • 000000 Li 

13 130.0000001 140. 000000 r 150.000000 r 

NPR = 5, 

APR (il_= - 

1 2. 000000 r 2.500000 r 3. 000000 » 

5 3.800000r 1.800000r 2.000000r 

9 2.4000007 2. 6000 00 > 2^800 0 

13 3.200000 r 3. 400000 r 3. 600000 r 

17 4.000000f . 4.200000r 4.400000r 


ARN (I)= 

1 i.000000f 1.000000r 0.700000f 

_ 5 0.^0 „j 0.0 f __ 

NSTG = 3r JCOOL = If 

DHQTD = 0.063500f ETATTD= 0.886000f 

TFFD = 51 . 529999 f XNRTD_= 40.100006f _ 

PSifi = 1.500000f XM2XD = 0.410000f 

ANGSWX= 2.899999r TTIN = 2500.00000 f 

PTIN „_r 59.800003f FARGD = 0. 02 0000 f 

END NAMELIST INPUT 


80. 000000 f 
80.000000f 
120.000000f 


3.500000 f 

2.200000f 

3. 000000 f 

^ . oeiavic>0 f 
4.600000f 


0.0 


NAMELIST INPUTl ^ _ 

KINDOF=f 864000f 

TC = 700.000000f FARCX = 0.0 f 

_YEAR = 1980.f.„ YEARB = 19^. f . 

ELIFE = 0.1000000E+05f 

END NAMELIST INPUTl 

PC;BLED= 0.06337 PC NCH= 0.020 0 1 E FF4= 0 .883 5 PRN= 3.315 „ 

NASA OUTPUT ON IFC=15fl6 

READY 


Figure 3. Input to Second Example Case (Concluded). 
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Table IV, Variable SeCtinga Co Trigger Default Calcu uCion 
of Some Deaig'.i Point Input. 


Variable 

Name 

Setting 

Action Taken 

NSTG 

0.0 

Program Calculates Number of Stages 
(Not Recommended) 

ETATTD 

0.0 

Program Calculates Design Point Efficiency 

PSID 

0.0 

Design Point Loading Set to 0.9 

ANGSWX 

180.0 

Exit Swirl Angle Calculated From Zero Hub 
React ion 

XMZXD 

0.0 

Sets Exit Axial Mach Number to 0.5 


Table V. Default Settings for Variables in Namelist "Input 1" 


Variable 

Name 



Units 

Default 

Value 

De.serlptlon 

KINDOF 

None 

86400000 

An Ordered Combination of Digits 
Representing the Cooling Configura- 
tion of the Turbine 

TC 

■’ R 

700.0 

Total Temperature of the Cooling 
Flow 

FARCX 

None 

0.0 

Fuel-Air Ratio of the Cooling Flow 

Year 

— 

1980 

First Year of Service for Stator 
Vane Material 

Year B 

— 

1980 

First Year of Service for Rotor 
Blade Material 

ELIFE 

hrs 

10,000 

Desired Life of Turbine Airfoil 
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The integer variable KINDOF represents the cooling configuration of the 
turbine. Each blade row starting with the first stage stator is assigned an 
integer value characterizing the type of cooling employed as follows: 

0 Uncooled 

1 Convection cooling 

2 Convection with coating 

3 Advanced convection 

4 Film with convection (755! trailing edge injection) 

5 Film with convection (505! trailing edge injection) 

6 Film with convection (25% trailing edge injection) 

7 Transpiration with convection (25% trailing edge injection) 

8 Full coverage film 

9 Transpiration 

For example, the 86400000 configuration has the first three blade rows 
cooled and the remaining five rows uncooled (a four-stage turbine). For a 
detailed description of the cooling flow calculation and the various cooling 
flow configurations, the reader should consult Reference 3. 
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4.0 PROGRAM OUTPUTS 


The basic output from the program consists of two tables. These tables 
show the turbine efficiency and turbine flow function variations for each 
of the first stage nozzle area ratios, pressure ratios, and pet: ant corrected 
speeds specified in the input. The input values of area ratio are converted 
to first stage nozzle angles before being printed out. The output tables for 
the first example case are shown on pap;.' , ’8 through 25. The table structure 
is compatible with NASA cycle deck reji'M ents given in Reference 2 (pages 
23 and 24). 


The output tables can be visualized as three dimensional, composed of a 
series of planes with each plane assigned a value of nozzle angle, BETA. 

Then in each BETA plane, the dependent variable (ordinate axis) is a function 
of pressure ratio, PR, and corrected speed, rpra. The dependent variables are 
respectively turbine corrected flow, W, and total-to-total efficienty, ETA. 


For example, in the output table on page 28, the six lines of the depen- 
dent variable correspond to the six values of corrected speed, respectively. 
And the five values of the dependent variable in each line corresponds to the 
five values of pressure ratio above each column, respectively. 


In addition to these two tables, there is a terminal listing summarizing 
the results of the cooling flow calculation, if this option was used. The 
value of the total cooling flow, PCBLED, is printed out together with the 
cooling flow for the first stage nozzle alone, PUNCH. The new cooled turbine 
efficiency value, EFF4, is given together with the new value of the total-to- 
total pressure ratio across the turbine, PRN. An example of this printout is 
shown on page 27 in the second example case. With the flows, shaft work, and 
turbine pressure ratio known, the user can calculate the new cooled turbine 
efficiency, ETATTD, using the bookkeeping procedure compatible with the cycle 
deck representation to be employed. A cycle deck efficiency scalar could 
then be used or, if desired, the program could be rerun on the uncooled branch 
using the new design point efficiency value as an input. 
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5.0 PROGRAM DIAGNOSTICS 


Hie PART computer program contains error printouts to aid the user in 
trouble shooting his input. A listing of the error messages and their mean- 
ings are given below. 

1. LIMITING VALUE OF UHUB-1 600.0, CALCULATED VALUE OF UHUB- 

This warning message is printed out only if the calculated rim speed 
exceeds the recommended value (this is a disk stress warning). 

2. LIMITING VALUE OF ANS-42.0E9, CALCULATED VALUE OF ANS- 

This warning message is printed out if the product of the exit annu- 
lus area and the rpm squared exceeds the recommended value (this is a cen- 
trifugal stress limit on the rotor blading). 

3. QIRE CTR ERROR— (CALLING LINE*, 15.) 

There are six iterations in the program. Each iteration is balanced 
using the Method of False Position. This method is contained in the sub- 
routine QIREXX. A maximum of 25 passes is allowed for any single iteration 
to balance. If the iteration does not balance within the specified toler- 
ance, the error message will appear with the number of the offending itera- 
tion in the 15 Format field. 

Normally, the occurrence of ,'>uch an error will not cause a problem. 
However, in the case of QIRE loop number five which calculates the turbine 
efficiency for the specified input values of pressure ratio and corrected 
speed, an additional message indicating the convergence error is printed out. 
This message has the form: 

DHQT* ,ERR* ,PQP- 

where the blanks contain the current values of specific enthalpy change 
divided by inlet total temperature, the convergence error in pressure ratio, 
and the pressure ratio at which the error occurred. 

The user should inspect the error to see if the degree of convergence 
is satisfactory, if not, it may be necessary to restrict the range of input 
speeds and/or pressure ratios requested. The individual QIRE loops together 
with the calling routines and type of iteration are as follows: 
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QIRE 

LOOP 

1 

2 

3 

4 

5 

6 


CALLING 

ROUTINES COMMENTS 

INLETX Calculates individual stage efficiencies from the input 
value of overall turbine efficiency (NSTG>1) 

VELRAT Obtains the axial velocity ratio across the rotor. 

CHOKEX Solves for the value of nozzle Mach number when the rotor 
chokes . 

ROTGKX Solves for exit annulus choke location given the location 
of rotor choke. 

PRTEFF Calculates efficiency for input values of speed and pressure 
ratio . 

FSTACK Solves for the "polytropic" exponent for a multistage tur- 
bine . 
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6.0 EXAMPLE CASES 


Two example cases are given in order to illustrate the use of the pro- 
gram. The first case utilizes the default settings to generate the output 
for a single-stage, uncooled, variable-geometry turbine. The second case is 
a four-stage, fixed-geometry turbine having cooling in the first three blade 
rows . 


A complete record of the two terminal sessions including a listing of 
the output tables is given on the following pages. The program inputs and 
outputs for these two cases have been discussed previously in Sections 3.0 
and 4.0. 


E14Z62A.FIl.E15. DATA DELETED 
E14262A.F1U16*DATA DELETED 
INPUT 

15r 


EXEC NPT 
ENTRY (A) 
.ENTRY lA) 
NAMELIST 
NSPDS * 
APCNC (I)« 
1 
5 

...... . 9 

13 

NPR « 

APR (D* 
1 
5 

.. 9. .. 

13 

17 


10.000000, 
50.000000, 
„90. 000000, 
130.000000, 
20 , 

1.100000, 

1.700000, 

2.400000.. 
3.200000, 

4.000000, 


20.000000, 

60.000000, 

100,000000, 

140.000000, 


1.200000, 

1.800000, 

2,600000, 

3.400000, 

4.200000, 


NAR 

ARN (I)« 
1 

NSTG » 
DHQTD s 
TFFD « 
PSID * 


ANGSMXs 
PTIN = 

END NAMELIST 


3, 


1 .30000b, 
ja*0 „.,i. 
1 , 

0.032780, 

62.979996, 

0.851100, 

15.200000, 

14.696000, 

INPUT 


1.000000, 

- . J).0 

JCOOL » 

ETATTD» 
XNRTD « 
XMZXD • 

TTIN ■ 
FARGO « 


30.000000, 

70.000000, 
110.000000, 
150.000000, 


1.400000, 
2.000O!?r.'< 
2 . 800000 , 
3.600000, 

4.400000, 


0.700000, 


0 , 

0.923000, 
193.520004, 
0.373000, 
518.669922, 
0.0 , 


40.000000, 

80.000000, 

120.000000, 


1.600000, 
2 . 2000(50 , 
3.000000, 
3.80000(5, 
4.600000, 


0.0 


ENTER NAMELIST INPUT ? 


& INPUT feEND 




NAMELIST 
NSPDS » 
APCNC (I)» 

INPUT 

15, 




1 

10 . 000000 , 

20.000000, 

30.000000, 

40.000000, 

5 

50.000000, 

60.000000, 

70.000000, 

80.000000, 

9 

90.000000, 

100.000000, 

110.000000, 

120.000000, 

13 

NPR * 

APR (1)» 

130.000000, 

20, 

140.000000, 

150.000000, 


1 

1.100000, 

1.200000, 

1.400000, 

1.600000, 

5 

1.700000, 

1.800000, 

2.000000, 

2 . 200000 , 

9 

2 . 400000 , 

2.600000, 

2:800000, 

3.000000, 

1.3 

_ ,3.200000, 

3.400000, 

3.600000, 

3.800000, 

17 

NAR 

ARN (I)= 

4.000000, 

3, 

4.200000, 

4.400000, 

4.600000, 

1 

5 

1.300000, 
0.0 , 

1.000000, 
0.0 , 

0.700000, 

0.0 , 

NSTG = 

1, 

JCOOL = 

0, 


DHQTD = 

0.032780, 

ETATTD= 

0.923000, 


TFFD = 

62.979996, 

XNRTD = 

193.520004, 


PSID = 

0.851100, 

XMZXD = 

0.373000, 


ANGSWX* 

15.199997, 

TTIN = 

518.669922, 


PTIN * 

END NAMELIST 

NASA OUTPUT 

14.696000, 

INPUT 

ON IFC=15,16 

FARGD = 

0.0 , 



READY 
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LIST FILE16.DATA 
FILE16.DATA 

270 TURBINE FLOW FUNCTION_VS._ PR» RPM»_AND BETA 

BETA 3 61 . 68. 75. 
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READY 
EXEC NPT 

£llZ^.EJLE15*nAlA. JELEIED 
ENTRY (A) E14262A.FILE16.DATA DELETED 


NAMELIST . 

INPUT 




NSPDS. « .. . 

15f . 

... ... 



APCNC (I)« 





1 

10.000000F 

20. 000000 F 

30. 000000 F 

40 

5 

_50,,000000i. 

60. 000000 F 

70. 000000 F 

80 

9 

90.000000F 

100. 000000 F 

110. 000000 F 

120 

r 13 . . 

130. 000000 F 

140. 000000 F 

150. 000000 F 


^ npr • ‘ _ 


... 5.3^ 



APR (T)* 





1 

1.100000F 

1. 200000 F 

1. 400000 F 

1 


_1, 7000001.. 

1. 800000 > 

2.000000 F 

2 

9 

2. 400000 F 

2. 600000 F 

2. 800000 F 

3 

13 

3.200000F 

3.400000F 

3. 600000 F 

3 

. _ 17 

§aJ000001__ 

4.200000F 

...4* 400000 F 

4 

NAR 

3 F 




ARN (D* 





1 

300000 r. 

1, 000000 F„ 

„„„ 0, 700000 F 

0 

5 

0.0 F 

0.0 r 



rJSTG ■ 

1 F 

irnni » 

^WWWW ” 

0F 


DHQTD « 

0.032780F 

ETATTD* 

0.923000F 


TFFD » 

62. 979996 F 

XNRTD * 

193.520004f 


PSID - 

0.851100F 

XMZXD « 

0.373000F 


ANGSWXs 

15.200000F 

TTIN * 

518.669922f 


PTIN « 

14.696000F 

FARGO » 

0.0 F 


END NAMELIST 

INPUT 





ENTER NAMELIST INPUT ? 

„ INPUT 

NSPDS«6 

APCNC«20 . » 40 . » 60 . » 80 . f 1 00 . » 1 20 . 

NPR^S 

APRa2.0»2.5»3.0r3.5»3.8 

NAR=1 

ARN=1.0 __ 

NSTG»3 

JCD0L»1 

DHQTD=0.0635 

ETATTD=.886 

TFFD=58.53 

XNRTD«40.10_ ■ 

PSID=1.5 
XMZXD*0.41 
.,„.ANGSWX=2 , 

TTIN=2500. 

PTIN«59.8' 

FARGDS0.02 _ _ 

&END 



000000 F 
000000 » 

000000 f 


600000 F 
200000 f 
000000 F 
800000 F 
600000 F 


.0 
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NAMELIST INPUTl 

J<INDOF»» 86400 (? 00 » ' 

TC* ‘ 700. 000000 » FARCX » " 0.0 

YEAR » . 1980.1 YEARB * 1980. r 

ELIFE * 0.1000000E+05r 

END NAMELIST INPUTl 


ENTER NAMELIST INPUTl J 
'■ INPUTl 
K I NDOF *864000 
SiEND 


MMILML_J,MPyi.„ 

NSPDS « 6» 

APCNC (I)« 


1 

20. 000000 r 

40.000000* 

60.000000* 

80.000000* 

5 

100. 000000 r 

120.000000* 

70.000000* 

80.000000* 

9 

90. 000000 » 

100.000000* 

110.000000* • 

120.000000* 

IS 

130. 000000 > 

140.000000* 

150.000000* 


NPR * 

APR n>* 

5» 




1 

2. 000000 r 

2.500000* 

3.000000* 

3.500000* 

5 

3.800000r 

1.800000. 

2.000000* 

2.Z00000* 

9 

Z.400000f 

2.600000* 

2.800000* 

3.000000* 

13 

3.200000* 

3.400000* 

3.600000* 

3.800000* 

17 

4.000000* . 

4.200000* 

4.400000* 

4.600000* 

, NAR * , . 

.1 f. 




ARN (I>= 





1 

1.000000* 

1.000000* 

0.700000* 

0.0 * 

5 

0.0 * 




NSTG « 

3* 

JCOOL = 

1 * 


DHQTD = 

0.063500* 

ETATTD= 

0.886000* 


TFFD = 

58.529999* 

XNRTD = 

40.100006* 


PSID = 

1.500000* 

XM2XD = 

0.410000* 


ANGSWX= 

2.899999* 

TTIN = 

2500.00000 * 


PTIN = 

59.800003* 

FARGD * 

__0, 020000* 


END NAMELIST 

INPUT 


• 


NAMELIST 

INPUTl 




KINDOF=f 

864000* 




TC 

700.000000* 

FARCX = 

0.0 * 


„JEAR ,= 

.1980. * 

YEARB = 1 780 . * 


ELIFE = 

0.1000000E+05* 




END NAMELIST 

INPUTl 





PQBLED = 0 . 063_37„,PXNCH^ 0 , t2JJi_lEF„4^__0 ,M35 . PRN = 1 , 31 5 

NASA OUTPUT ON IFC=15rl6 

READY 
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LIST FILE15.DATA 
FILE15.DATA 

Z69 TURBINE EFFICIENCY^AZS. PRr RPM» AND BETA 
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LIST FILE16.DATA 
F1LE16.DATA 

270 TURBINE FLOW FUNC TION VS. PR» RPM» AND BETA 

BETA 1 71. 






• 1 
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7.0 ANALYTICAL BACKGROUND 


The following section has been written in order to give the user a gen- 
eral idea of the type of turbine representation used in the program and the 
approach used in the derivation of the equations. Details of the deriva- 
tions together with sample calculations may be found in the Monthly Progress 
Reports (e.g., References 10, 11, and 12). 


7.1 TURBINE MAP REPRESENTATION 


Typically, cycle deck entry to a turbine map is through corrected speed, 
N/SQRT(T), and actual energy, DH/T, with turbine flow function, W*SQRT(T)/P, 
and total-to-total efficiency being output. Total-to-total pressure ratio is 
sometimes used instead of actual energy as the second map entry. 

The discussion of the turbine map representation can be conveniently 
subdivided into two parts: the flow and the efficiency. 


The flow model is illustrated by the three sketches shown in Figure 4. 
The two curves on the top in the figure are used to generate the flow repre- 
sentation on the bottom. Sketch 4-1 shows the turbine stage characteristic 
(i.e., a plot of turbine loading against flow coefficient). Sketch 4-2 shows 
the dependence of the maximum value of the turbine flow function on corrected 
speed. With the corrected speed and DH/T known, the stage loading can be 
calculated, and the flow coefficient obtained from Sketch 4-1. Once the flow 
is choked (i.e., the choked branch of the stage characteristic), the flow 
coefficient remains constant for that speed. Sketch 4-2 is next used to ob- 
tain the maximum value of the turbine flow function at the corrected speed of 
interest. The value of the turbine flow function is then calculated. The 
equations used are as follows: 


U/At - (21TR/60) (N//rRggoT 
- (DH/T)/[2(U//T)2/goJ] 
C/A{. ■ ((})/Cos a 2 )/(U/A(.) 


TFF 


(TFF) 


max 
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where 

DH/T >• turbine stage specific enthalpy drop divided by inlet total 

temperature, Btu/lbra R 

N/SQRT(T) « cor ected speed, RPM/SQRT(R) 

R * pitch line radius, ft 

C * velocity, ft/ sec 

Aj. ■ speed of sound at inlet total temperature, ft/sec 

IFF ■ inlet turbine flow function, lbm/sec*SQRT( R)/psia 

Cos(a 2 ) ■ nozzle exit angle 

The value of Co8(o2) is obtained from the design point information as 
is the value of the pitch line radius. 

The efficiency model is illustrated by the five sketches shewn in Figure 
5. The first three curves In the figure are used to generate the last two 
sketches. The "backbone" of the turbine map shown in Sketch 5-4 is the locus 
of the peak efficiency at each value of DH/T. This locus is obtained from 
^ Sketch 5-1, which shows turbine pitch line loading along the map "backbone" 

as a function of DH/T. The "backbone" efficiencies are obtained from Sketch 
I 5-2. This sketch gives the "backbone" loss (defined as the difference be- 

tween the ideal and actual values of DH/T) as a funciton of DH/T. With the 
"backbone" loading and efficiencies known at each DH/T, Sketch 5-3 is used to 
I evaluate the "off-backbone" loss and to obtain the efficiency at any value 

of corrected speed. Wlien the turbine "off-backbone" loss is plotted with the 
coordinates shown in Sketch 5-3, the resulting curves are nearly linear at 
any given DH/T. These five curves, three univariate and two bivariate, are 
r sufficient to define the turbine map. Note that as shown on Sketch 5-5, the 

design point does not generally fall on the "backbone" but is separated from 
it by a "stand-off" distance which is calculated from design point informa- 
tion. 

I 

I The analytical basis for the five correlating curves is discussed in the 

following sections. 

7.2 TURBINE FLOW MODEL 

f The turbine stage characteristic serves as che basis for modeling the 

turbine flow. An analytical expression for the stage characteristic of a 
; constant-pitch, axial-flow turbine can be obtained by using the continuity, 

energy, and angular momentum equations, together with a number of relation- 
ships from the pitch line vector diagram. In deriving this equation, it is 
assumed that the pitch line flow angles at nozzle and bucket exit are in- 
I variant, and that the axial velocity ratio across the bucket is constant. 


32 








This lacer asaumpcion ia reaaonable for incompreaeible or low Mach number 
(M<0.25) flowB. It ahould be emphasized, however, that the final correla- 
tions have no such restrictions. 


The stage characteristic for a single-stage turbine with the above as- 
sumptions can be written in the form 

ip " i^(tan 02 tan 03)- 1/2 (1) 

where 

t(i “ DH/(2U^/goJ) , Turbine stage Loading 
4 i ■ Cg/U, Flow coefficient 

02 “ Nozzle exit flow angle 

03 ■ Bucket exit flow angle. 

Selecting a reference point and assuming that 


tan 02 


^z3 a 

+ tan 63 

Cz2 


a constant, 


this equation may be written as 

JL_ = i_ + _1_ (2) 

tR 4>R 2% \4 ir / 

For unchoked flow, the reference point (indicated by the subscript R) was 
selected at the design point on the turbine map. For choked flow, the refer- 
ence point was selected at the critical point (subscript GRIT) on the turbine 
map. The critical point is located at the value of DH/T at the reference 
speed (e.g., 100%) where nozzle choking first occurs. For a variable-geometry 
turbine, the angle 03 is calculated from the input value of nozzle area ratio. 

The equation used for the normalized flow coefficient is 

^ _ C^/At (N//t)r _ C/At (N//Tr 

h " (Cz/At)j^ N//f (C/At)j, N//r 

The velocity ratio, C/Aj-, for a point at a selected speed is determined 
by calculating a pseudoarea at which the Mach number is assumed to be unity 
(i.e., at the maximum turbine flow function for that speed). This pseudoarea 
is assumed to be constant for that speed. Tlie Mach number (velocity ratio) 
at any point on the speed line is then calculated from the usual flow function 
equations. By definition, the pseudoarea varies with speed in direct propor- 
tion to the maximum turbine flow function variation with speed. This method 
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of calculating the flow coefficient wae found to give better reaulte than 
that obtained using the first stage nozzle throat area for all speeds. For 
the case in which the first stage nozzle is choked, the two procedures are 
identical . 

Typically, the calculated and measured values of velocity ratio are 
with in about 6%, with the larger error#, occurring at the higher Mach numbers . 
The predicted value of the normalized flow coefficient intercept (at DH/T*0) 
are within about 5% of experimental ly derived valuits . These errors combine 
to yield flow errors on the order of 5 % for nominal, area maps ( i.e. , nozzle 
area ratios equal to 1). Slightly higher values may occur for other stator 
settings. 

A typical comparison between measured and calculated values of the tur- 
bine flow function is shown in Figure 6. This figure is for test turbine 
number 25 in Table I. The maximum error shown on this figure is about 2.4% 
and occurs at the lower end of the test data. This type of plot was used in 
obtaining the error estimates given above. 


7.3 TURBINE LOSS MODEL 


the 


There are four key steps in the devslopmsnt of 
turbine off-design loss model. These steps are 


Cu6 €v{Uatt.ons governing 


1. The development of an equation giving the turbine total-to-total 
efficiency at a general point in terras of nozzle and bucket effi- 
ciencies coupled with a seraiempirical loss term due to the depar- 
ture of the rotor incidence angle from the optimum. 

2. The transformation of Che seraiempirical incidence angle loss law 
so as to eliminate the explicit occurrence of the incidence angle 
by introducing the stage loading. 

3. The differentiation of the resulting efficiency expression in 
order to obtain the locus of peak efficiencies. This peak effi- 
ciency ridge Chen becomes the "backbone" of the map. 

4. The substitution of the peak efficiency relationships back into 
Che general efficiency equation in order to obtain an expression 
for Che ’’off-backbone" loss. 


The development of the efficiency expression proceeds from the (h,s) 
diagram for adiabatic flow through a two-dimensional turbine stage as shown 
in Che following sketch. 
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Comparison of Measured and Calculated V.alues of Inlet 
Turbine Flow Function. 
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Enthalpy-Entropy Diagram for a Turbine Stage 


Usir.g the station numbers shown In the above sketch, the definition of 
turbine total-to-tota] efficiency may be written In the form: 


rvp'j' 


AH 


AH + 


To 3 

T 3 


(h3 - h3s) 


+ 


^o3s 

T2 


(h2 - h2g) 


(A) 


By definition, the nozzle efficiency is equal to the ratio of the actual 
nozzle exit kinetic energy to the ideal nozzle exit kinetic energy. A similar 
definition in terms of relative velocities holds dor the bucket efficiency. 

For off-design calculations at any incidence angle (i), an additional loss 
term must be included, i.e., 


h3 . -coshi))^ 

where 

Htq = nozzle efficiency 

Hg *' bucket efficiency 

W2 “ bucket inlet relative velocity 

W3 = bucket exit relative velocity 
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The semiempirical incidence angle loss law is based on the assumption 
that the kinetic energy of the component of velocity normal to the optimum 
incidence angle is lost. This is a fairly standard assumption (see, for 
example, Reference 4). Powers other than 2 are frequently used on the 
cosine . 


Before introducing the incidence angle loss term into the efficiency 
exnression, it was transformed into the following expression 



The stage characteristic was then used to substitute for flow coefficient 
in terras of stage loading. The substitution of this results into the expres- 
sion for efficiency yielded, after simplifying, the following results. 


where 


1 

ViTT 

- 1^ ■ A (2ij) + 

1)2 

+ B (ijj - 




(7) 


J3 Uq V ' 



To3g 

/i_ - 1 

) ^ 


, i 

V^z2, 

/ Cos2s3 

T2 

V TIN 

/ Cos2a2_ 

(8) 

’ 4 


tan 

z '■'z2 

tan 33 

2 



'^03 S2op 

Tq 2 (2^op + 1)2 


The temperature ratios in Equations 8 and 9 are of order one as i' the 
bucket axial velocity ratio. If 'hese ratios together with the blad • iw 
efficiencies and exit flow angles are assumed to remain constant then Equa- 
tion 7 can be differentiated, and the peak efficiency point located. 



( 10 ) 


__1 ^ ^ A (2 + 1)^ -f B 

^pk 'I'pk 


Note that if the nozzle and bucket efficiencies equal unity in Equation 
(9), then A = 0.0. '"hen there is no loss other than incidence and = iji^p 
as expected. 
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By substituting Equations 10 and 11 into Equation 7, the following ex- 
pression for "off-backbone" loss can be obtained. 

Equations 10, 11, and 12 give the location of the peak efficiency, the 
magnitude of the peak efficiency, and the variation in efficiency as we move 
away from the peak. 

These equations represent the stage loss characteristic of a turbine. 
The design point information is used to obtain the initial values of A and B 
as well as the values of the blade row efficiencies and metal angles. The 
loss equations are then applied at incremental values of DHQT starting at 
zero to obtain the turbine efficiencies. Approximate relationships are used 
for temperature and velocity ratios to obtain new values of A and B for each 
DHQT. 


Typically, the calculated values of the loss slopes and those obtained 
from air turbine test data are within about 5% for corrected speeds within 
plus or minus 20 % of the design point value. The values of the "backbone 
efficiencies" generated by the above equations do not include either Reynold' 
Number effects or Che severe rotor exit losses encountered near exit annulus 
choke. In order to account for these effects, the loss along the peak effi-’ 
ciency ridge was empirically modified using the results of the NASA air 
turbine tests. Although relatively good correlation existed between the 
different test turbines at low values of DHQT, the drop in efficiency in the 
neighborhood of exit annulus choke was so severe that correlation was diffi- 
cult. For this reason, variations in efficiency on the order of 5 % can be 
obtained in this region, 

A comparison between measured and calculated values of the turbine total 
to-total efficiency is shown in Figure 7. Tliis figure is for test turbine 
number 25 in Table I. The maximum error shown on this figure is about 0.8%. 
This type of plot was used in obtaining the error estimates given above. 
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Turbine Total-to-Total Efficiency 



Figure 7. Comparison of Measured and Calculated Values of Total-to-Total 
Efficiency. 
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LIST OF SYMBOLS 


A, a 

B. b 
C 

Cp 

go 

h 

1 
J 
N 
P 

Rg 

s 

T 

U 

W 

a 

B 

AH 

r 

X 

P 

0 

'P 

<P 

n 

Subscripts 

1 

2 
3 
B 
d 
N 
0 

op 

pk 

TT 

S 

z 


Constants 
Constants 
Velocity, fps 

Constant Pressure Specific Heat, Btu/lbm*R 
Dimensional Constant, 32.17 ft Ibin/lbf sec^ 

Enthalpy, Btu/lbm 

Incidence Angle (i * 62 “ ^2op^ » degrees 
Mechanical Equivalent of Heat, 778.16 ft. Ibf/Btu 
Speed, rpm 
Pressure, psia 

Gas Constant, 53. 35 ft lbf/lbm“R 

Entropy, BtuAbm’R 
Absolute Temperature, ’ R 
Wheel Speed, fps 
Relative Velocity, fps 

Angle of Absolute Velocity With Axial, degrees 

Angle of Relative Velocity Vector With Axial, degrees 

Drop in Total Enthalpy, Btu/lbm 

Ratio of Specific Heats 

Enthalpy Loss Coefficient 

Fluid Density, Ibm/ft^ 

Ratio or Total Temperature to Standard Temperature 
Turbine Stage Loading [t = AH/ ( 2U^/gQj) ] 

Flow Coefficient (iji = C 2 /U) 

Ef f ic iency 


Nozzle Inlet 
Rotor Inlet 
Rotor Exit 
Bucket 

Design Point Value 

Nozzle 

St agnat ion 

Opt imum 

Peak 

Total-to-Total 

Isentropic 

Axial 
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